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I.  Introduction 


An  understanding  of  the  fundamental  Interactions  between  a  solid  surface 
and  adsorbed  atoms  or  molecules  is  necessary  for  a  proper  description  of  many 
Important  processes  such  as  catalytic  reactions,  corrosion,  vapor  deposition 
and  gas-surface  scattering.  Theoretical  studies  along  these  lines  range  from 
classical  and  semlclasslcal  to  purely  quantum  mechanical.  While  most  have 
addressed  the  somewhat  more  tractable  problem  of  weak  pbysl sorption, 1  some 
models  for  strong  chemisorption  have  also  appeared.  In  the  present  work  we 
are  concerned  with  the  vibrational  states  of  the  system,  namely  the 
vibrational  structure  of  the  adspecles,  nf  the  bond  between  the  surface  and 
the  adspecles,  and  that  of  the  bulk  solid  and  Its  surface  (the  phonons).  In 
the  past,  the  dynamical  Interactions  between  the  solid  and  the  adspecles  have 
Eluded  multiphonon  effects,  one-phonon  effects  or  no  phonons.  There  Is 
strong  Indication  7rom  some  of  these  studies  that  multiphonon  effects  are 
quite  Important,  a  particular  example  being  stimulated  desorption.3  The  use 
of  high-power  lasers  to  Initiate,  control  or  enhance  gas-phase  reactions  is 
established.4  Lasers  have  also  been  used  as  heat  sources  on  surfaces 

in  applications  such  as  annealing5  and  deposition  via  thermal  reactions  in 
the  gas  phase.6  Only  a  few  applications  of  low-power  lasers  to  gas-surface 
reactions  have  been  reported,  however,  and  selective,  nonthermal  laser 
stimulated  surface  processes  (LSSP)  seem  to  be  taking  place  In  some  of  these 
experiments.7'8  The  present  work  Is  an  attempt  to  set  up  a  general 
formalism  capable  of  dealing  with  LSSP  from  a  microscopic  point  of  view  which 
Is  as  "first  principles"  as  practicable. 

Recent  calculations  of  gas-surface  Interactions  within  a  scattering  or  de¬ 
sorption  context  have  adopted  classical  and  semlclasslcal9  or  quantum3,10,11,12 
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approaches  with  varying  degrees  of  phenomenology.  The  semi classical  approaches, 
for  example,  represent  a  good  compromise  between  the  simplicity  of  classical 
trajectories  and  the  details  available  from  full  quantum  treatments.  The  latter 
tend  to  be  computationally  tedious,  often  provide  less  Insight  and  require  approxi¬ 
mations  such  as  simplified  phonon  band  structure,  one-phonon  treatment  with  a 
linear  chain  model  of  the  solid,  and  neglect  of  memory  effects  in  the  statistical 
reduction  of  the  lattice  many-body  problem. 

The  gas-surface  studies  referred  to  In  the  previous  paragraph  have  dealt 

with  either  physlsorption  or  chemisorption,  but  no  single  formulation  seems  to 

be  applicable  v.o  both.  In  the  phenomenological  theories  of  chemisorption,  the 

strong  cheml sorptive  bond  cannot  be  treated  perturbationally.  A  unitary  trans- 

13  14 

formation  of  the  unperturbed  states  Is  necessary  to  account  for  this.  *  The 
states  of  the  adatom  are  thereby  changed  to  those  of  quasi  particles  which  move 
on  the  surface,  carry  a  distortlonal  wave  along  with  them,  and  Include  most  of 
the  atom-surface  Interaction  energy. 

Another  approach  Is  to  Ignore  the  two-dimensional  periodic  nature  of  the 
surface  and  Include  only  a  few  surface  atoms  close  to  the  adatom.  Such  "cluster" 
calculations  give  rise  to  discrete  energy  levels  in  contrast  to  the  band  type 
energy  level  structures  typical  of  systems  with  t-anslational  symmetries.  The 
significance  of  the  band  structure  In  gas -surface  dynamics  Is  one  of  the  ques¬ 
tions  which  we  address  here. 

The  gas-surface  Interaction  Is  a  complex  many-body  problem.  In  treatments 

of  the  dynamics  of  surface  processes,  classical  and  quantum  statistical  theories 
3  0-13  15-18 

have  been  advanced  ’  *  to  explain  the  general  qualitative  nature  of  the 

thermodynamics  (desorption,  migration,  etc.).  However,  quantitative  Interpreta¬ 
tion  Is  difficult  because  fundamental  parameters  such  as  the  energy  relaxation 
rate  and  the  surface  migration  rate  are  not  easily  available  from  experiment; 
they  can  be  obtained  theoretically,  however,  the  accuracy  being  limited  by  the 
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type  of  potential  used  (Morse,  LEPS,  Lennard- Jones ) ,  the  degree  of  dimensionality, 
and  simplifying  assumptions  such  as  additivity  of  transition  rates  due  to  different 
forces.  In  the  case  of  LSSP,  we  further  need  to  Include  the  (possibly 
synergistic)  effects  of  the  absorption  of  coherent  electromagnetic  radiation 
by  the  adsorbed  species. 

A  few  experimental  studies  of  LSSP  have  been  reported.  7» 8  Large 
enhancements  of  the  rates  of  reactions  (relative  to  the  rates  for  thermally 
Initiated  reactions)  Involving  hydroxyl  and  amino  groups  on  a  silica  surface 
have  been  obtained  using  low-power  (-10  W/em?)  C02  laser  radiation. 7  These 
reactions  take  place  on  the  surface,  but  the  exact  mechanism  responsible  for 
the  enhancement  is  unclear.  Efficient  selective  laser  pumping  of  particular 
vibrational  inodes  of  the  adsorbed  species,  followed  by  dissociation  and 
desorption  of  OH  and  NH  groups  have  been  suggested.  Alternatively,  a 
mechanism  Involving  a  "hopping"  electron-transfer  mechanism  has  been  proposed. 

While  both  mechanisms  seen  plausible,  there  Is  no  conclusive  evidence  from 
theoretical  or  experimental  studies  to  support  either.  An  estimated  very 
short  lifetime  (-  10*11  s  )  0f  the  vlbratlonally  excited  states  of,  say,  the 

7 

S1-N  bond  would  tend  to  favor  the  electron- transfer  mechanism.  However, 
lifetimes  as  short  as  these  do  not  appear  universally,  and  one  cannot  dismiss 
the  vibrational  pumping  mechanism.  Other  Instances  of  LSSP  have  been 
reported,  some  Involving  absorption  of  the  radiation  in  the  gas  phase  with 
subsequent  surface  reaction,  and  others  Involving  absorption  of  radiation 

a 

during  an  actual  surface  process. 

Recent  theoretical  attempts  aimed  at  understanding  the  mechanisms  In  LSSP 
have  adopted  a  mainly  statistical  kinetic  approach.  3,7»17»18  They  have 
pointed  out  the  sensitivity  of  LSSP  to  the  laser  frequency  and  fluence,  and 
have  shown  that  the  critical  parameter  Is  the  rate  of  energy  transfer  from  the 


vlhratlorally  excited  adspecles  to  the  solid.  This  energy  relaxation  rate 

In  turn  depends  on  the  nature  of  the  coupling  between  the  adspecles 

solid.  While  lattice  vibrations  represent  the  primary  means  for 

energy  transfer  In  nonmetals.  In  metals  the  quasicontinuum  of  conduction 

electrons  can  exchange  energy  with  excited  vibrational  states  of  the 
19 

adspecles.  Evaluation  of  these  rates  seems  to  be  the  primary  bottleneck  In  a 
quantitative  understanding  of  LSSP. 

In  this  paper,  a  formalism  Is  presented  to  deal  with  the  energy 
relaxation  of  a  vlbratlonally  excited  adatom  pimped  by  an  Infrared  laser  and 
coupled  to  phonons  of  the  solid.  This  represents  only  a  small  portion  of  the 
vast  array  of  processes  oceurlng  In  a  laser-irradiated  gas-surface  system,  hut 
extension  to  other  aspects  of  the  problem  Is  possible. 

The  core  of  this  formalism  Is  presented  In  Section  II.  The  discussion  In 
Section  III  Includes  limiting  cases  and  suggests  points  of  departure  where 
contact  can  be  made  with  other  formalisms,  mainly  phenomenological  ones.  The  main 
points  of  the  paper  are  simmarlzed  In  Section  IV. 


II.  Formalism  for  Multiphonon  Relaxation  of  an  Adatom  After  IR  Laser 

Irradiation. 

The  physical  processes  which  this  formalism  attempts  to  describe  Involve 
the  absorption  of  Infrared  laser  radiation  by  an  adatom  bound  by  the  periodic 
(In  2D)  potential  at  a  solid  surface.  The  vlbratlonally  excited  state  thus 
created  then  decays  via  energy  transfer  to  the  various  modes  of  vibration  of 
the  solid.  For  the  sake  of  simplicity  we  do  not  consider  relaxation  to 
electronic  degrees  of  freedom  which  would  be  quite  important  for  metals. 
Furthermore,  it  Is  assumed  that  a  good  representation  of  the  potential  at  the 
surface  can  be  obtained  either  In  coordinate  space  or  In  the  Fourier  space. 

For  example.  In  many  cases,  particularly  physlsorption,  the  surface  potential 
can  be  approximated  by  a  sum  of  pair  potentials, 10  the  Fourier  expansion 
of  which  converges  rapidly  enough  to  make  the  problem  computationally 
tractable.  The  fact  that  only  energy  transfer  to  phonon  modes  of  the  solid  Is 
being  considered  means  that  adatom-adatom  Interactions,  which  lead  to  Tz'  or 
phase- relaxation  processes,  are  being  Ignored.  In  other  words,  a 
low-coverage  sltutatlon  Is  being  envisaged. 

In  considering  the  Interaction  between  the  laser  and  the  adatom,  previous 
models  have  often  assumed  the  existence  of  a  well-defined  "active" 
mode,  17,18,20 that  Is  one  associated  with  a  transition  between  states  of  the 
adatom-solid  system  which  are  separated  by  an  energy  gap  nearly  equal  to  the 
laser  frequency.  In  the  Bloch  picture  of  delocalized  adatom  states,  the 
energy  levels  appear  as  bands,  with  band  gaps  which  are  complicated  functions 
of  the  20  wave  vectors  In  the  2D  Brill ou in  zone  of  the  surface.  The 
Identification  of  a  particular  single  active  mode  is  difficult  In  this  ease. 

In  dealing  with  direct  transitions  of  the  adatom  from  states  In  one  band  to 
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those  In  another,  the  "resonance"  condition  may  occur  at  many  different  parts 
♦ 

In  the  20  k-space,  and  It  1$  the  cumulative  effect  of  all  such  transitions 
which  wUT  distinguish  one  laser  frequency  from  another.  For  the  moment  the 
only  restriction  that  may  be  applied  to  the  laser  Is  that  Its  frequency  be  of 
the  same  order  of  magnitude  as  the  expected  gaps  In  the  energy  level  band 
structure  of  the  adatom  under  the  Influence  of  the  surface  potential. 

The  other  Interaction  being  considered  Is  that  between  the  vlbratlonally 
excited  adatom  and  the  phonon  modes  of  the  solid.  If  adsorption  Is  not  very 
strong,  one  may  proceed  by  solving  the  zero-order  problem  of  the  adatom  In  t 
potential  field  of  the  solid  frozen  In  Its  equilibrium  configuration,  and 
applying  the  Influence  of  the  motion  of  the  lattice  particles  on  the  adatom  as 
a  perturbation.  With  zero-order  states  of  the  system  represented  by 
combinations  of  the  above  zero-order  states  of  the  adatom  and  phonon  states  of 
the  harmonic  lattice  In  the  absence  of  the  adatom,  multiphonon  relaxation  of 
the  adatom  Is  described  by  transitions  between  these  maqy-quaslparticle 
states.  The  strength  of  these  transitions  Involves  matrix  elements  of  the 
perturbation  operator  mentioned  before. 

A.  The  Total  Hamiltonian. 

*#0 

Let  J(a  be  the  Hamiltonian  of  the  adatom  of  mass  m  under  the  Influence 
of  the  potential  V°(r)  due  to  the  solid  at  equilibrium, 

Kfc).-—  V*  ♦  V°(?).  (1) 

2m  r 

V°(r)  -  IV(r-R°).  (2) 

1  * 

7 
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♦ 

where  r  Is  the  position  vector  of  the  adatom  and  R“  Is  the  equilibrium 
position  of  the  1th  lattice  particle.  Here  V(r-R^)  Is  a  suitably  chosen  pair 
potential  between  the  adatom  and  the  X  lattice  particle.  The  actual 
potential  can  be  written  as 

V(r)  -  IV(r-R).  (3) 

l  * 

■*  th 

where  R^ls  the  Instantaneous  position  of  the  l  lattice  particle.  The 

♦  -  «♦ 

difference  between  V(r)  and  Vu(r)  Is  the  perturbation  tM.  responsible  for 
phonon  relaxation  of  the  vlbratlonally  excited  adatom, 

AJC  -  V(r)  -  V0(r)  -  l[v(r-R  )  -  V(r-R°)]  •  (4) 

1  ^  X 

To  get  the  conplete  Hamiltonian  for  the  system,  the  contributions^  andj(r 
due  to  vibrations  of  a  harmonic  lattice  and  interaction  with  the  laser  field, 
respectively,  must  be  Included*  *•«•» 

x  -ks  *xr  ♦  .  (M 


He  employ  solutions  of  the  two  zero-order  problems,  namely  one  set  for  the 
adatom  In  the  potential  vO(r)# 


(6) 


and  another  set  for  the  normal  modes  of  vibration  of  -,he  lattice  In  the 
harmonic  approximation  , 


with 


uO  0  0  0 

y  Yy  *  Cy  Yy 


(7) 


xj  ■  (1/2)  l  {qf(k)  ♦  .  (8) 

♦  D  DO 

kb 
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♦ 

In  Eqs.  (6)-(8),  a,  j  and  n  are  quantum  numbers  representing  a  band  Index, 

quantization  of  motion  along  the  z  axis  and  quantization  of  motion  In  the  xy- 

♦  0  ♦ 
plane,  respectively  (n  Is  a  2D  wave  vector);  the  eigenvalues  correspond 

n  0  0 

to  the  stationary  states  of  the  adatom  In  the  potential  V°;  yv  ®nd  *v  are  the 
many-phonon  wave  function  and  energy,  respectively,  for  the  whole  crystal;  and 

4  *  -4 

q  (k)  and  q  (k)  are  the  normal  displacement  coordinate  and  momentun, 
b  b 

respectively,  for  a  phonon  with  wave  vector  k  In  the  phonon  branch  b  and 

frequency  u  (k). 
b 

The  Hamiltonian  for  the  Interaction  of  the  system  with  the  laser 
field  Is  simplified  by  assuming  that,  due  to  a  near-resonance  condition, 
most  of  the  radiation  Is  absorbed  by  the  surface-adatom  bond.  The  Hamiltonian 
Xr  then  becomes 

Hr  m  K  (P*A  *  A*p),  (9) 

2m  0 

«♦ 

where  p  Is  the  momentum  operator  for  the  adatom,  e^  is  some  residual  effective 
charge  on  the  adatom  and  A  Is  a  vector  potential  for  the  laser  field  at  time  t, 

A  (r,t  )  *  Ao  ?  exp  [1(ic.r  -  ut  )]  +  c.c.  ,  (10) 

c,  k  and  u  being  the  unit  polarization  vector,  wave  vector  and  frequency, 
respectively  .  If  it  Is  further  assumed  that  the  Incoming  flux  of  photons 
Is  large  enough  so  that  absorption  dominates  over  emission,  then  the  complex 
conjugate  part  (c.c.)  of  Eq.  (10)  can  be  neglected. 


B.  Solutions  of  the  Zero-Order  Problem. 

The  zero-order  problem  Is  represented  by  the  Hamiltonian 

jc°  -ks  o» 
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For  simplicity  the  assumption  thatJCj  andjfy  can  be  solved  separately  [Eqs. 
(6)  and  (7)]  is  made.  Hie  solutions  so  obtained  are  then  combined  to  provide 
a  description  of  the  overall  problem  defined  ^  the  Hamiltonian  In  Eq.(5). 

The  two  zero-order  problems  are  discussed  below. 


..0 

1.  Zero-Order  Problem  for  the  Adatom. a,. 

As  In  most  problems  Involving  periodicity,  It  Is  convenient  to 

transform  to  the  Fourier  space.  At  the  surface  of  a  perfect  solid,  the 

potential  Is  periodic  along  the  surface  plane  (thexy-plane)  but  not  In  the  z 

0  - 

direction.  The  appropriate  Fourier  expansion  of  V  (r)  Is  based  on  the 


condition 


0  ♦  *0  0  * 
v  (r  ♦  X^)  -  V  (r). 


where  Is  a  lattice  vector  along  the  surface  plane 

♦0  0  0  ♦  ♦ 

X  n,a.  *  n„a„, 

l  ll  2l  11  2  2 


02) 

03) 


•»  ■* 

ni,  n2  being  Integers  and  ai,  82  primitive  surface  lattice  vectors.  As  a 

0  ♦ 

result  of  Eq.  (12)  we  can  write  V  (r)  as  a  Fourier  series 

0  ♦  0  ♦  0  .!  ♦ 

V  (r)  s  V  (X.z)  •  l  (z)  e1G-x  .  (14) 

6 

0 

where  the  Fourier  transform  V*  (z)  Is  given  by 

G  r 

ve <z) "  4;  x  dx  (i5) 

% 

2D 

unit  cell 

♦  ♦ 

Au  being  the  area  of  a  20  unit  cell.  Ay  ■  |a-|  x  821*  The  expansion  In  Eq.  (14) 
Is  over  all  20  reciprocal  lattice  vectors (RLV)  6  for  the  exposed  surface  plane, 

10 


2  ■  *  ®2^2 


b,  - 


2  IT  (»2  X  2) 

<s2  *  *) 


2  t 


2v  (2  x  a^) 


(16) 

(1?) 


<a2  x  2) 


“1  *  '“2  "  *'  1 
where  2  is  a  unit  vector  normal  to  the  surface  (xy-  or  aja2-) 
plane.  Thus  for  any  plane  parallel  to  the  surface  plane,  the  2D 
RLV  are  completely  specified  by  two  integers  and  a2  to  be  used 
as  subscripts,  for  example  Vg  <*> 

The  surface  potential  Is,  In  general,  a  smoothly  varying  function, 

and  the  Fourier  expansion  can  be  expected  to  converge  quite  rapidly.  This  Is 

borne  out  by  recent  calculations#1*3,21  using  summed  pair  potentials  (Morse, 

Yukawa  etc.),  where  (2)  was  found  to  be  appreciable  only  for  mlt  mg  <2, 

and  by  the  qualitative  success  of  models  using  pure  sinusoidal  functions  for 

V°(r).  22  The  maximum  value  of  (mi .mg)  can  therefore  serve  as  one  measure  of 

the  degree  of  approximation  In  actual  computations. 

The  stationary  states  for  a  single  particle  with  a  periodic  Hamiltonian 

I/*3  are  characterized  by  a  quasi  continuum  of  quantixn  numbers,  namely  20  wave 
Jla 

vectors 


♦  V1  ♦  vz  * 

h  *T  bl 


(18) 


"1  "2 

v1  and  v2  are  Integers,  and  N]  and  Ng  signify  the  (large-)  number  of  unit 

cells  along  the  ai  and  a2  directions,  respectively,  beyond  which  the  crystal 

0 

surface  Is  repeated.  The  single-particle  states  e  "kist  then  satisfy  the 
Bom-von  Kaman  boundary  condition  on  a  surface 


0  ♦  *  0  ■*  *  0  ♦. 

*  (r  ♦  Niai)  ■  *  (r  ♦  ^2)  •  *  (r). 


(19) 


leading  to  the  Bloch  functions 


0  ♦  .  0  ,  .  16-X  1n*X 

♦.j;  (r)  ‘  l  *  • 

G 
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(20) 


In  Eq.{20),  summing  over  G  allows  us  to  restrict  n  ■  ♦  n2  to  the  2D 

first  Brlllouln  rone  (FBZ), namely  the  area  defined  by 


bi  bi  b2  b2 

~  <  Di  <  tl2  <  — 

Uslng  Eqs.  (1),  (14)  and  (20)  In  Eq.  (6),  we  get 


-  f  f>2  2  .  0  1G'*X  0  ( 

♦  I  2m  r  ♦  G'  «Jn 

G  L  G'  J 


(21) 


X  exp  (1G*X)  exp  (1q*X)  ■  0.  (22) 
Using  symmetry  associated  with  translation  of  the  reciprocal  lattice  by  RLVt 
we  obtain 


_(T  1!?  2  A2  >  *  2  0  1  0 

jlralVS  <’'*G)  *  ‘.jsJ 


(*> 


0  0 

I  V+  ♦  (*)  ♦ 

♦  t  G-G’  •jn-Kr 


w 

► 

« 


(r)  [•  exp  (1(n  ♦  G)*X}  •  0 


(23) 


which  leads 
namely. 


to  a  set  of  coupled  differential  equations  for  the  e  ♦ 


[H2  <j2  *2  *  *  •  0  J  0 

2m  ^z2  *  2m  ^  *ajnj 


(*) 


(*). 


*  l  yLe- 

G* 


(r)  ■  0. 


(24) 


The  solution  of  the  set  of  equations  (24)  represents  the  first  major 
computational  task  In  a  treatment  of  the  stationary  states  of  an  adsorbed 
atom. 


12 


4/0 

2.  Zero-Order  Problem  for  the  Lattice  VI brat Ions. u. 


h 


The  transformation  of  the  actual  Hamiltonian  of  a  vibrating 

lattice  to  the  form  of  Eq.  (8)  Is  based  on  a  Born-Oppenhelmer  (adiabatic) 

approximation,  and  the  Inclusion  of  only  harmonic  terms  In  the  expansion  of 

the  Interatomic  potential  a  ;  here  £  denotes  a  particular  electronic  state  of 

s 

the  crystal.  The  normal-mode  transformation  reduces  the  vibrating  lattice  to 
a  set  of  3pN  noninteracting  one-dimensional  harmonic  oscillators,  where  a  Is 
the  number  of  atoms  In  a  unit  cell  and  N  Is  the  number  of  unit  cells  In  the 
crystal. 

In  principle,  we  can  write  Eq.  (7)  as  a  Schrodinger  equation, 

* 


ii 

kb 


-  V 


pq2  + 

%  (k) 


2  *  2  * 

♦  Vk)Vk)J 


v  ■ 


0  0 
l  7 
V  V 


(27) 


and  write  as  a  product  of  single-oscillator  wave  functions 


l<k)]  ■ 
•CD 


(28) 


si  nee  X  is  a  sum  of  Independent  single- oscillator  Hamiltonians,  with  vt>(k) 
v  * 

the  number  of  quanta  of  the  (kb)  vibration  In  the  system.  Eq.  (27)  can  now  be 
written  as  a  series  of  single-oscillator  aquations 

0 


4 


H*  b2  12  ♦  2 

-7-77  *T“»  '"’"b  00 

®Mk) 

D 


vi(k) 


tv‘>] 


■  «°  *  ♦  fq  (k)l  , 

vb(*)  VO 1  b  J 

whose  solutions  ara  £v  ■  v^k),  q  *  q^k),  w  «  «b(k)J 


(29) 
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*°(q) 

v 


-(1/2)«q2/h 

v!  e 


H 

v 


where  H 

v 


Is  the  v 


th 


Hermlte  polynomial. 


Also, 


0 

c 

V 


(30) 


(31) 


where  the  subscript  {v}  represents  a  specific  set  (v  (k)}  of  3pN  vibrational 

b 

quantum  numbers  which  describes  the  state  of  the  system. 

The  dispersion  relation  embodied  In  the  functional  form  of 

«  (k)  Is  the  primary  Information  necessary  for  a  complete  description 
b 

of  the  phonon  spectrum.  Extensive  experimental  and  theoretical  work  has  been 
done  in  this  area.and  highly  accurate  phonon  band  structures  and  dispersion 
relations  are  known  from  neutron  scattering.  Infrared  absorption,  Raman 
scattering,  Brlllouln  scattering  and  X-ray  scattering.  23 


3* 

In  the  harmonic  approximation,  the  equilibrium  positions  of  the  atoms  In  a 
*0 

crystal,  R^,  do  not  represent  dynamical  variables.  The  total  zero-order 
Hamiltonian,  Eq.  (11),  can  then  be  written  as 


1  2*  Z.  , 

2  ub(k)S(lc) 


h  . 


(32) 
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The  solutions  of  a  Hamiltonian  which  Is  the  sum  of  Independent  single-particle 
Hamiltonians,  In  this  case  for  one  adatom  and  3pN  phonons,  are  products  of 
solutions  of  the  Individual  single-particle  Schrodlnger  equations. 


0 


a°.«.  (r)"["[  4»°  ♦  [q  (k)"]. 

•to  ♦  %(k)  L  b  J 

DK 


(33) 


These  will  be  denoted  |  ajri;  (v  (k)}  >,  and  the  (r,q)  representation  will  be 

b 

Implied.  Also,  the  abbreviated  forms  v,  q  and  u  will  be  used  Instead  of 
v  (k),  and  «h(k),  respectively,  wherever  appropriate. 

C.  Absorption  of  Laser  Radiation  and  Phonon  Relaxation. 

Eq.  (5)  Is  the  total  Hamiltonian  for  the  system.  One  possible  way  of 

viewing  It  Is  to  ascribe  to  &H  and  Hr  the  status  of  perturbations  Imposed  on  a 

system  otherwise  described  by  the  zero- order  Hamiltonian}!  »K#  +X^  •  Such 

a  description  Is  not  appropriate  for  all  cases  of  physical  Interest,  and  it  Is 

necessary  to  clearly  define  the  limitations  of  such  an  approach  as  well  as 

alternate  routes  for  cases  falling  outside  these  limits. 
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The  recent  phenomenological  studies  of  LSSP  *  have  provided  some 
Insight  Into  the  above  questions.  It  Is  seen  that  relative  time  scales  for 
the  vibration  of  the  adatom/surface  bond  and  the  phonons,  the  associated 
energy  gap  and  the  effective  coupling  between  the  two  must  be  considered  In 
deciding  what  kind  of  perturbative  approach.  If  any.  Is  applicable.  Thus,  for 
high  temperatures  (large  bX),  high  optical  pimtplng  rates  (large Kr),  very 
strong  bonding  between  the  adatom  end  the  surface,  or  a  very  small  energy 
difference  between  the  adatom/surface  vibrational  frequencies  end  the  Debye 
frequency  of  the  solid,  straightforward  application  of  a  perturbation  approach 
may  be  Inadequate.  In  such  cases  it  may  become  meaningless  to  separate  the 


i 

* 

4 

\ 

i 


J 


degrees  of  freedom  of  the  solid  from  those  of  the  edetom/surface  bond,  and  a 
basis-set  transformation  may  be  necessary  to  describe  the  new  degrees  of 
freedom  which  would  In  general  differ  considerably  from  the  original  ones. 

Such  transformations  have  been  discussed  previously  in  the  context  of 

phenomenological  treatments  within  a  second-quantized  formulation  of 
phonon-mediated  energy  transfer  In  condensed  systems.  We  shall  not  pursue 
this  aspect  of  the  problem  but  restrict  ourselves  to  cases  where  &Xandj(r 

4r0  4.O 

can.  In  fact,  be  considered  as  perturbations  of 

We  begin  by  writing  the  transition  rate  for  an  adatom  from  a  laser-pmped 
state  |e>  of  energy  cc  to  a  state  |g>  of  energy  eg  using  Fermi's  Golden  rule 

>  -Bcu-1  -Mi  -Be, 

*(9*e)  •  (T«  )  (•  «e  -  •  (34> 

X  |<1e|  AX|fg>|2  6  [ec  -  eg  -  (e*-ei)]  , 
wherell)andl1>are  Initial  and  final  phonon  states  of  the  lattice, 
respectively,  with  energies  «i  and  c^,  end  W9  and  Wg  are  probabilities  of 
findings  the  adatom  In  the  states  |e  >  end  |g  >,  after  a  fast  short  pulse  of 
laser  radiation.  Because  of  the  6-functlon,  we  must  have 

•  *i  •  cp  •  eg  •  Aeg  "  ^“eg  •  (35) 

which  allows  us  to  write  It  as 

X|<1e|  AJl|fg>|2  lUeg-  (e#  -  ei)]  f  (36) 


or 


IKg*e)  ^  JJ  e  ^weg  |<le|  aJtlfgsf  *  [A*g-(*j-eiJ.  (37) 


To  simplify  H(9*-e)  further  It  Is  convenient  to  expend  attln  terms  of  Fourier 
components  of  the  pair  potential  v(r).  With 


_  ♦  * 
f  ♦  -1q*V 

V**T.  V(r)«  &  • 

q  \) 


where  V  is  the  volume  of  a  unit  cell#  we  can  write 
u  ♦  <*  *0 
0  ♦  1q*(r-R  ) 

V  (r)  -  l  1 

*  Q 
l  R 


hr) -II . 
♦  q 
i  q 


so  that 


r  w'r  I  -Ir«»  . 

1^3  *  I  V  "  J 


The  rate  then  becomes 


*'’•*> •  £ ($* ' )  H*' 


•*  ♦  ♦,  ♦  ♦  ♦  ♦  j0 

1q*r  *  -1q'*r  -1q*R*  -1q’** 

<#lv^  •  «  !•>  •  J  !,<*!•  -  •  K> 


♦  ♦  0  1 
1q'**i*  Iq'**!'  | 

X  <f|e  ®e  |l>U(atfl  -  C«f  ®*i3) 
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We  now  use  the  6 -function  representation 


In  the  operator  which  correspond  to  the  product  of  two  exponentials  will 
survive.  The  remaining  terms  Involve  operators  ©  (other  than  unity)  which 
create  or  annihilate  at  least  one  phonon  and  contribute  nothing  to  <1|6M>* 
The  constant  terms  Involved  In  the  phonon  operator  lead  to  6(oegt/h)  when  the 
Fourier  transform  Is  taken,  and  these  also  do  not  contribute.  The  overall 
Fourier  transform  can  therefore  be  written  as 


m 

I 


dt  e 


lAggt/h 


«< 


-Iq.u  <t) 
e  * 


1q-u  .(•) 


»• 


(49) 


The  double  brackets  represent  an  ensemble  average  whose  evaluation  has  been 

24  + 

discussed  extensively  In  the  literature.  The  operator  u^  Is  linear  In  the 
creation  and  annihilation  operators  for  phonons.  Hence  It  follows  that  the 
comnutator  [A.B3  Is  a  c-nunber,  where  A  and  8  are  the  exponents  In  the 
ensemble  average,  that  Is 


A  1  q*u  (t)  , 

(50) 

B  1  }.  0^,(0)  # 

(51) 

Thus 

*B  .  *A*B  t  (1/2)[A.B] 

(52) 

and 

« 

(t)  -  q'*u  ,  (0)  ] 

»  •  «e  1  1  ». 

(53) 

For  operators  b  linear  In  harmonlc-osclllator  creatlon/annlhllatlon 
25 


operators  we  have 


1#  (-1/2)«6Z»_ 
«e  »  •  o  T* 


(54) 


Use  of  the  above  eventually  leads  to  the  result 


11,26,27 


J 

a 


! 

4 

•? 

i 
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where 


« 


»  m  .XP^,  (q.q1)  ♦  q’C^.ft)  *  qj. 


*•£•^*9')  s  (1/2)  <<  (^.u. )  ^  +  (q1 


U‘  •  -£l 

Is  a  Debye-Waller-1 Ike  factor,  and 


(55) 


(56) 


CiA,(t)  s  «UA(t)uA. (0)» 


(57) 


Is  a  correlation  function  Involving  the  atomic  displacements  at  sites  1,1* 
and  times  t  and  0.  The  energy  transfer  rate  thus  becomes: 


“  1  HVjVj.n  U,  * eg 

q  q  q»eg  <1  ,eg 

<1  <1 


0  0 

-1(q.R  -q'.R  ) 

•  II,  •  *  1 

i  l* 


-«o 


The  separation  of  the  rate  expression  Into  sections  depending  on  the  adatom 

properties  and  those  depending  on  the  properties  of  the  solid  Is  an  attractive 

feature  of  this  formalism  and  is  primarily  due  to  the  treatment  of  momentum 

28 

transfer  and  energy  transfer  as  Independent  variables.  The  computational 
procedure  now  calls  for  the  solution  of  the  coupled  differential  equations  for 
the  stationary  states  of  the  adatom  In  the  static  crystal  potential, 
calculation  of  transition  probabilities  for  these  states  under  laser  ^ 
Irradiation,  and  finally  the  calculation  of  the  correlation  functions  C  (t). 
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One  starts  by  utilizing  the  expansion  of  the  displacement  u^, 

in  the  harmonic  case,  in  terms  of  phonon  creation  and  annihilation 
t  29 

operators  a£b  and  a^,  respectively.  With  M  the  mass  of  the 

lattice  atoms  in  a  monatomic  crystal  and  N  »  the  total 

number  of  unit  cells,  we  can  write 

to  obtain  the  correlation  tensor  In  the  form 

Cv'1  v'0^'  $?rJ 

kb  «k(k) 

— 

1 k.(R  -R  )[  -V(k)t  _  1u(k)t 

*  e  til  (7^^  +  1)  e  b  +  n^  e  b  ,  (60) 


where  n*.  is  the  Bose  function 


"*b“ 


(  pfMk)  Vl 

K  b  7  • 


It  Is  expedient  at  this  point  to  take  the  continuum  limit  of  the  phonon 

♦ 

band  structure  fay  replacing  the  stnmatlon  over  wave  vectors  k  with  an 


Integration  over  k  space. 


where  V  is  the  volume  of  the  crystal,  to  obtain 
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C„.(t)  *  J E£u_ I  dk 

“  T&S  b 


4  4  ♦  r  ♦  "1 

„  1k*(R  -R  J  -V(k)t  ,  ♦  I 

(k)eh*(k)  e  *  *je  *>  +2n^bCos{ub(k)tJ  .1 


From  Eq.  (58),  the  rateft(9+e)  Is  seen  to  Involve  the  time  Fourier  transform 

of  exp[C  ,(t)].  If  a  cutnul ants-type  expansion  Is  reasonable,  the  first-order 

**  29 

term  will  be  Important,  namely 


|  dt  e1Ue9t  CiA,(t)  -  |  n(*eg)  +  '}  sgn  (u>gg) 

I  ^  r*  •*  ♦  *  +  1k.(R  -R  .)  2  2  I 

x  I  I  dk  j  eb(k)eb  (k)  e  1  x  M“eg  -  wb(k)}J  * 


The  term  Involving  Cw,(t)  Is  the  n-phonon  term.  Clearly  the  one-phonon  term 


(n-1)  dominates  when  «  1,  where  a  Is  a  typical  lattice 

parameter.  Evaluation  of  the  higher-order  terms  can  be  simplified  by  writing 


the  rate  as  a  Taylor  expansion  In  the  C^tt).  We  first  wr1teH(g*e)  as 


4U9*e) 


•(S’vfc 


1uftgt  ^ 

dt  e  9  I  2  (C  ,(t))  . 

IX 'W  XX 


where  *  Is  the  characteristic  time  of  the  highest  frequency  phonon.  Writing 
0 
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‘6CW*$5'' 


(66) 


(t,5')  being  the  Cartesian  Indices  x,y,z,  we  get  the  rate  In  the  form' 

I  l  i_  •  o  (h»g )/_  a"?  \  , 

U  J  n»l  JU'  Iff  n»*t  I  I 

Y Cti' *55'/  0 


(67) 


where 


0n,U'.W'(“*sl 


1 


luegt  n  ,  . 
dt  e  C  ,  ,  t)  . 

JU'  KK 


(68) 


Use  of  the  relationship 


(^)F«-Wl  (*2)  -  (^1^2)  '  . 


(69) 


where  Fw(f)  Is  the  Fourier  transform  of  f(t)  at  u. 


F“(,)  ‘  k 


idtf(t)e 


lut 


(70) 


gives  a  recursion  relation  between  the  *fs: 
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0* 

(»e3>  •  L  \  ...  »!  (7,) 

l,u',tC  ]  "•**  •'  ’  '  - 


which  is  useful  for  computational  purposes. 
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III.  Discussion 

In  developing  Eq.  (58)  for  the  vibrational  energy  transfer  rate  from  a 
laser-excited  adatom  to  phonons,  the  two-dimensional  periodicity  along  the  xy- 
plane  has  not  been  explicitly  exploited.  Important  physical  effects  and 
simplifications  of  the  formalism  are  associated  with  this  periodicity,  some  of 
which  we  now  examine. 

0  ♦  ■* 

The  Potential.  The  nature  of  the  potentials  V  (r)  and  V(r)  In  Eqs.  (2), 
(3),  (39)  and  (40)  Is  very  different  along  the  z  direction  as  compared  to  that 
along  thexy-plane.  Both  must  approach  a  constant  for  large  z  and  be 
essentially  oscillatory  on  the  xy-plane  near  the  surface.  It  may  also  be 
necessary  to  allow  for  different  contributions  to  the  potentials  due  to  Ions 
at  or  near  the  surface  and  those  deep  In  the  bulk,  formally,  the  latter  Is 
accomplished  by  affixing  an  extra  Index  to  the  pair  potential  In  Eqs.  (2)  and 
(3),  namely, 

vV)  •  I  *  (72) 

Jl  *  *  1 

V(r)  •  f  r-R .)  (73) 

A  *v  *  • 


Inclusion  of  these  effects  Is  probably  most  practical  for  purposes  of  a 
phenomenological  description  of  LSSP.  Writing  Eqs.  (39)  and  (40)  In  the  forms 


and 


+  *0 
1Q»X 

A 


(74) 


♦o  ♦  ♦  ♦  ♦  ♦ 

the  sunmatlon  over  X^  In  Eq.  (74)  reduces  to  6(Q-G),  with  q  *  (Q,qz)  end  G  a 

primitive  2D  reciprocal  lattice  vector  In  the  plane  of  the  surface.  Such  a 

♦ 

reduction  does  not  take  place  In  Eq.  (75)  because  the  X  are  not  fixed 

i 

vectors.  Because  of  the  smallness  of  the  atomic  displacements  from 

equilibrium,  It  Is  reasonable  to  use  a  sharply  peaked  (Gaussian  or  lorentzlan) 

function  fc*(Q) ,  with  the  peak  at  Q  ■  G,  to  replace  the  summation  over  X  in 
G  Jt 

Eq.  (75).  Eqs.  (74)  and  (75)  may  therefore  be  written  as 


V°(r) 


l  l 

q  ♦ 
Q 


and 


. . 

6(Q-G)e  *  1 


♦  ♦ 
1q*r 


(76) 


V(r)-  I  l  l 
q  j  z 
z  Q  l 


♦  -1q  z 

*.(Q)e  *  * 

G 


1q*r 


(77) 


The  term^g(Q)  represents  the  "Inelasticity"  of  the  transition  process  viewed 
as  a  scattering  event,  namely  the  departure  from  purely  diffractive  scattering 

4  4 

described  by  6(Q-G)  In  t<\.  (76).  The  averaging  procedure  Implied  In  the  use 

of&  +  (Q)  Is  not  always  appropriate.  In  particular,  if  the  characteristic  time 
G 

scale  of  the  overall  relaxation  process  Is  conparable  to  or  smaller  than  a 
typical  period  of  vibration  of  the  atoms  In  the  solid,  no  such  averaging  Is 
possible.  Examination  of  Eqs.  (42)  -  (46)  with  the  aim  of  applying  the  forms 
(76)  and  (77)  of  the  potential  reveals  a  substantial  simplification  by 
allowing  a  reduction  of  the  sunmatlons  over  R^  and  R^,  to  those  over  only  z^ 
and  z^,.  This  Is  an  Important  step  In  establishing  the  connection  between 
three-dimensional  and  one-dimensional  treatments  of  ISSP.  The  development 
of  a  formalism  to  exploit  this  two-dimensional  translational  symmetry  Is  In 
progress. 
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-1/2  *  1k*R 

Polarization  Eigenvectors.  The  terms  N  f(k)e  1 

O 


In  Eq.  (59)  represent  components  of  the  eigenvectors  of  the  dynamical  matrix 


In  the  harmonic  approximation. 


lhe  z-component  of  k  must  clearly  be 


complex  In  order  for  the  displacement  function  to  damp  out  for  z»0.  The 
displacement  operator  may  then  be  written  as 


-  /*y/2r „  „  t(k  £)***■  iK.xr  ti 

"*'W  1 1  **  [  *  (y‘-T 


Assuming  a  knowledge  of  the  form  of  u  (k),  the  expression  In  curly  brackets 
♦  b 

depends  only  on  b,  K  and  z^,  all  dependence  on  kz  having  been  absorbed  In  the 

summation.  This  represents  another  point  at  which  a  phenomenological  approach 

+ 

could  be  applied  effectively,  since  a  model  for  u  (k)  could  be  Introduced 

b 

along  with  assumptions  regarding  the  range  of  values  for  kz. 

♦  ♦  4 

Models  for  Dispersion.  In  addition  to  e.(k),  a  knowledge  of  u  (k)  Is 

b  b 

necessary  for  actual  evaluation  ofA(g»e).  Detailed  Information  on  the  form 
♦ 

of  u  (k)  Is  available  from  experimental  and  theoretical  studies  of  a  number  of 
b 

systems,  but  for  the  present  semi quantitative  study,  a  model  such  as  the 
Einstein  or  Debye  model  may  suffice.  In  the  former,  phonon  dispersion  Is 
entirely  suppressed  by  assuming  the  form 

u  (k)  ■  u  ,  a  constant,  (79) 

>  b  D 

which  neglects  all  correlations  between  displacements  of  different  lattice 

♦ 

atoms,  end  ty  ascribing  a  single  vibrational  frequency  to  all  points  k  on  all 
branches  b.  This  model  precludes  energy  transfer  over  a  range  of  energies 
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but  does  provide  the  simplest  picture  of  phonon  band  structure.  The  Debye 
model,  on  the  other  hand,  assures  a  linear  dispersion 

»b(k)  •  k  (80) 

Independently  of  branch  Index,  and  thereby  allows  for  energy  transfer  over  a 

(quasi-)  continuous  range.  The  use  of  Eq.  (80)  In  conjunction  with 

assumptions  of  Isotropy  and  simple  cubic  structure  reduces  the  expression  for 

n 

C^,(t)  and  the  Fourier  transforms  of  [C^( (t)]  to  computationally  convenient 
forms.  Actual  calculatlonal  procedures  are  currently  being  developed  In  this 
laboratory. 


ZB 


IV.  Summary  and  Conclusions 

A  general  theoretical  technique  has  been  developed  for  the  treatment  of 
wultlphonon  relaxation  of  an  adatom  vlbratlonally  excited  by  low-potter  laser 
radiation.  The  dynamical  processes  are  embodied  In  correlations  between 
displacements  of  different  atoms  of  the  solid  at  different  times.  The 
correlation  functions  and  their  Fourier  transforms,  though  complicated, 
provide  a  concise  physical  description  of  the  relaxation  process;  the 
n-phonon  processes  involve  Fourier  transforms  of  the  nth  power  of  the 
correlation  function.  The  generality  of  the  final  expressions  for  the 
relaxation  rate  has  the  advantage  of  providing  links  between  phenomenological 
treatments  of  the  problem  and  establishing  the  significance  of  various 
assumptions  Involved  In  the  latter. 


,9 


Acknowledgements 


We  would  like  to  thank  Paul  DeVries,  William  Murphy  and 
Michael  Hutchinson  for  helpful  discussions.  Special  thanks  are 
due  Horia  Metiu  for  a  critical  and  constructive  reading  of  the 
manuscript.  This  work  was  supported  in  part  by  the  Office  of  Naval 
Research  and  the  Air  Force  Office  of  Scientific  Research  (AFSC) , 

United  States  Air  Force,  under  Grant  AFOSR-82-0046.  The  United  States 
Government  is  authorized  to  reproduce  and  distribute  reprints  for 
governmental  purposes  notwithstanding  any  copyright  notation  hereon. 
TFG  acknowledges  the  Camille  ar>d  Henry  Dreyfus  Foundation  for  a 
Teacher-Scholar  Award  (1975-82) . 


30 


References 


[1]  W.  A.  Steele,  Surface  Scl.  36,  317  (1973);  F.  0.  Goodman,  CRC  Crlt.  Rev. 

Solid  State  Mater.  Scl.  7,  33  (1977). 
m  G.  Doyen  and  G.  Ertl,  J.  Chem.  Phys.  68,  5417  (1978). 

[3]  J.  I.  Gersten,  R.  Janow  and  N.  Tzoar,  Phys.  Rev.  B  11_,  1267  (1975); 

H.  J.  Kreuzer  and  0.  N.  Lowy,  Chein.  Phys.  Lett.  78,  50  (1981);  0.  Lucas  and 
G.  E.  Ewing,  Chem.  Phys.  58,  385  (1981);  C.  Jedrzejek,  K.  F.  Freed,  S. 

Efrlma  and  H.  Metlu,  Surface  Scl.  109,  191  (1981);  G.  Korzenlewski ,  E.  Hood 
and  H.  Metlu,  J.  Vac.  Scl.  Technol.  20,  594  (1982). 

[4]  T.  F.  George,  I.  H.  Zlmnerman,  P.  L.  DeVries,  J.-M.  Yuan,  K.-S.  Lam,  J.  C. 
Bellum,  H.-W.  Lee,  M.  S.  Slutsky  and  J.  Lin,  In  Chemical  and  Biochemical 
Applications  of  Lasers,  vol.  IV,  edited  by  C.  B.  Moore  (Academic  Press, 

New  York,  1979);  articles  on  Laser  Chemistry  by  A.  H.  Zewail,  V.  S.  Letokhov, 
R.  N.  Zare,  R.  B.  Bernstein,  Y.  T.  Lee  and  Y.  R.  Shen  In  Physics  Today  33, 
(11),  25-59  (1980). 

[5]  G.  K.  Celler,  J.  M.  Poate  and  L.  C.  Kimerllng,  Appl.  Phys.  Lett.  32,  464 
(1978);  A.  B.  Fowler  and  R.  T.  Hodgson,  Appl.  Phys.  Lett.  36,  914  (1980); 

A.  Bhattacharyya  and  B.  G.  Streetman,  J.  Phys.  D  L67  (1981). 

[6]  C.  P.  Christensen  and  K.  M.  Lakln,  Appl.  Phys.  Lett.  32,  254  (1978). 

[7]  M.  S.  Djidjoev,  R.  V.  Khokhlov,  A.  V.  Kiselev,  V.  I.  Lygin,  V.  A.  Namiot, 

A.  I.  Osipov,  V.  I.  Panchenko  and  B.  I.  Provotorov,  In  Tunable  Lasers  and 
Applications,  edited  by  A.  Mooradlan,  T.  Jaeger  and  P.  Stokseth  (Springer, 
Berlin, 1976);  M.  S.  Dzhldzhoev,  A.  I.  Osipov,  V.  Ya.  Panchenko,  V.  T. 
Platonenko,  R.  V.  Khokhlov  and  K.  V.  Shaitan,  Zh.  Eksp.  Teor.  Flz.  74,  1307 
(1978)  [English  translation  In  Sov.  Phys.  JETP  47,  684  (1978)]. 

[8]  S.  D.  Allen  and  M.  Bass.  J.  Vac.  Scl.  Technol.  16.,  431  (1979);  T.  F.  Deutsch, 
D.  J.  Ehrlich  and  R.  M.  Osgood,  Jr.,  Appl.  Phys.  Lett.  35,  175  (1979); 

D.  J.  Ehrlich,  R.  M.  Osgood,  Jr.,  and  T.  F.  Deutsch,  Appl.  Phys.  Lett.  36, 

916  (1980). 

31 


[9]  S.  A.  Adelman  and  J.  D.  Doll,  J.  Chem.  Phys.  61_,  4242  (1974);  J.  D.  Doll, 

L.  E.  Myers  and  S.  A.  Adelman,  J.  Chem.  Phys.  4908  (1975);  Y.-W.  Lin 
and  S.  A.  Adelman,  J.  Chem.  Phys.  68,  3158  (1978);  M.  Shugard,  J.  C.  Tully 
and  A.  Nitzan,  J.  Chem.  Phys.  6£,  2534  (1977);  J.  C.  Tully,  J.  Chem.  Phys. 
73,  1975  (1980);  R.  I.  Masel ,  R.  P.  Merrill  and  W.  H.  Miller,  J.  Chem. 

Phys.  64,  45  (1976). 

[10]  F.  0.  Goodman  and  W.-K.  Tan,  J.  Chem.  Phys.  5£,  1805  (1973);  F.  0.  Goodman 

and  I.  Romero,  J.  Chem.  Phys.  6£,  1086  (1978);  G.  Wolken,  Jr.,  J.  Chem. 

Phys.  6fl,  2210  (1974);  Y.-W.  Lin  and  G.  Wolken,  Jr.,  J.  Chem.  Phys.  65, 

2634  (1976);  F.  J.  Milford  and  A.  D.  Novaco,  Phys.  Rev.  A  4,  1136  (1971); 

H.  Chow  and  E.  D.  Thompson,  Surface  Scl.  59^  225  (1976);  W.  E.  Carlos, 

G.  Derry  and  D.  R.  Frankl ,  Phys.  Rev.  B  1£,  3258  (1979);  W.  E.  Carlos  and 

M.  W.  Cole,  Phys.  Rev.  Lett.  43,  697  (1979);  S.  Mukamel ,  J.  Chem.  Phys. 

70,  2479  (1979). 

[11]  H.  Metlu,  J.  Chem.  Phys.  67,  5456  (1977). 

[12]  S.  Efrlma  and  H.  Metlu,  J.  Chem.  Phys.  69,  2286  (1978);  S.  Efrima,  K.  F. 
Freed,  C.  Jedrzejek  and  H.  Metlu,  Chem.  Phys.  Lett.  74,  43  (1980);  G.  Doyen, 
Phys.  Rev.  B  22,  497  (1980);  see  also  Reference  [3). 

[13]  K.  Kitahara,  H.  Metlu,  J.  Ross  and  R.  Sllbey.  J.  Chem.  Phys.  65^,  2871 
(1976). 

[14]  M.  Grover  and  R.  Sllbey,  J.  Chem.  Phys.  5£,  4843  (1971). 

[15]  J.  L.  Beeby  and  L.  Dobrzynskl,  J.  Phys.  C  4,  1269  (1971). 

[16]  H.  J.  Kreuzer,  Surface  Scl.  100,  178  (1980). 

[17]  T.  F.  George,  J.  Lin,  K.-S.  Lam  and  C.  Chang,  Opt.  Eng.  19,  100  (1980); 

J.  Lin  and  T.  F.  George,  Surface  Scl.  100,  381  (1980). 


32 


[18]  J.  Lin,  A.  C.  Berl,  M.  Hutchinson,  W.  C.  Murphy  and  T.  F.  George, 

Phys.  Lett.  79 A,  233  (1980);  T.  F.  George,  A.  C.  Berl,  K.-S.  Lam  and 

J.  Lin,  In  Laser  Applications,  vol.  5,  edited  byR.  K.  Erf  and  J.  f. 
Ready  (Academic  Press,  New  York),  in  press. 

[19]  J.  W.  Gadzuk  and  H.  Metlu,  Phys.  Rev.  B  22,  2603  (1980);  G.  Korzenlewski , 

T.  Maniu  and  H.  Metlu,  Chem.  Phys.  Lett.  73,  212  (1980);  J.  Chem.  Phys. 

76,  1564  (1982);  H.  Metlu  and  J.  W.  Gadzuk,  J.  Chem.  Phys.  2i»  2641 
(1981);  S.  Efrlma  and  H.  Metlu,  Surface  Sci.  109,  109  (1981);  J.  W.  Gadzuk 
and  H.  Metlu,  In  Vibrations  at  Surfaces,  edited  by  R.  Caudano,  J.  M.  Giles 
and  A.  A.  Lucas  (Plenum  Press,  1982),  p.  519  ff. 

[20]  A.  Nitzan  and  J.  Jortner,  Molec.  Phys.  25,  713  (1973). 

[21]  A.  Tsuchlda,  Surface  Sci.  14,  375  (1969);  see  also  Reference  [10]. 

[22]  R.  I.  Masel,  R.  P.  Merrill  and  W.  H.  Miller,  J.  Chem.  Phys.  65,  2690  (1976); 

N.  Garcia  and  N.  Cabrera,  Phys.  Rev.  8  18,  576  (1978);  N.  Garcia,  F.  0. 

Goodman,  V.  Celll  and  N.  R.  Hill,  Phys.  Rev.  B  19,  1808  (1979). 

[23]  P.  A.  Egel staff,  editor.  Thermal  Neutron  Scattering  (Academic  Press,  New 
York,  1965);  G,  Venkataraman,  L.  A.  Feldkamp  and  V.  C.  Sahni ,  "Dynamics 
of  Perfect  Crystals"  (MIT  Press,  Cambridge,  Massachusetts,  1975),  Ch.  5 
and  references  therein. 

[24]  C.  B.  Duke  and  G.  E.  Laramore,  Phys.  Rev.  B  2,  4765  (1970). 

[25]  A.  Messiah,  "Quantum  Mechanics,"  vol.  I  (John  Wiley  &  Sons,  New  York, 

1958),  Ch.  XII. 

[26]  C.  Klttel,  "Quantum  Theory  of  Solids"  (John  Wiley  &  Sons,  New  York,  1963), 
Ch.  19. 

[27]  B.  Bendow  and  S.-C.  Ying,  Phys.  Rev.  B  7_,  622  (1973). 

[28]  L.  Van  Hove,  Phys.  Rev.  95_,  249  (1954). 

[29]  A.  A.  Maradudln,  E.  W.  Montroll,  G.  H.  Weiss  and  I.  P.  Ipatova,  "Theory 

of  Lattice  Dynamics  In  the  Harmonic  Approximation,"  2nd  edition  (Academic 

Press,  New  York,  1971),  p.  63, 


33 


SP472-3/A1 


472: GAN:  715:enj 
78u472-608 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  GZN 


No. 

Copies 

Office  of  Naval  Research 

Attn:  Code  472 

800  North  Quincy  Street 

Arlington,  Virginia  22217  2 

ONR  Western  Regional  Office 

Attn:  Or.  R.  J.  Marcus 

1030  East  Green  Street 

Pasadena,  California  91106  l 

ONR  Eastern  Regional  Office 
Attn:  Dr.  L.  H.  Peebles 
3uildlng  114,  Section  0 
666  Sumner  Street 

Boston,  Massachusetts  02210  1 

Director,  Naval  Research  Laboratory 
Attn:  Code  6100 

Washington,  D.C.  20390  1 

The  Assistant  Secretary 
of  the  Navy  (RE&S) 

Department  of  the  Navy 
Room  4E736,  Pentagon 

Washington,  D.C.  20350  l  • 

Commander,  Naval  Air  Systems  Command 
Attn:  Code  3  ICC  (H.  Rosenvasser) 

Department  of  the  Navy 

Washington,  D.C.  20360  1 

Defense  Technical  Information  Center 
3uilding  5,  Cameron  Station 
Alexandria,  Virginia  22314  12 

Dr.  Fred  Saalf eld 
Chemistry  Division,  Code  6100 
Naval  Rastarch  Laboratory 
Washington,  D.C.  20373 

Or.  David  L.  Nelson 

Chemistry  Program 
Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  Virginia  22217 


No. 

Copies 


U.S.  Army  Research  Office 
Attn:  CRD-AA-I? 

P.0.  Box  12211 

Research  Triangle  Park,  N. C,  27709  1 

Naval  Ocean  Systems  Center 

Attn:  Mr.  Joe  McCartney 

San  Diego,  California  92152  1 

Naval  Weapons  Center 
Attn:  Dr.  A.  3.  Ams  tar , 

Chemistry  Division 

China  Lake,  California  93555  l 

Naval  Civil  Engineering  Laboratory 

Attn:  Dr.  R.  W.  Drisko 

Port  Rueneme,  California  93401  1 

Department  of  Physics  A  Chemistry 

Naval  Postgraduate  School 

Monterey,  California  93940  1 


Scientific  Advisor 
Commandant  of  the  Marine  Corps 
(Code  RD-1) 

Washington,  D.C.  20380 

Xeval  Shi?  Research  and  Development 
Center 

Attn:  Dr.  G.  Bosaajian,  Applied 
Chemistry  Division 
Annapolis,  Maryland  21401 

Navel  Ocean  Systems  Center 
Attn:  Dr.  S.  Yamamoto,  Marine 
Sciences  Division 
San  Diego,  California  91232 

Mr.  John  3oyie 

Materials  3ranch 

Naval  Ship  Engineering  Center 

Philadelphia,  Pennsylvania  19112 


1 


SP472-3/A17 


472: GAN: 716. -lab 
78u472-608 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  056 


No. 

Copies 

Dr.  G.  A.  Somorjai 

Department  of  Chemistry 

University  of  California 

Berkeley,  California  94720  1 

Dr.  L.  N.  Jarvis 
Surface  Chemistry  Division 
4555  Overlook  Avenue,  S.W. 

Washington,  D.C.  20375  1 

Dr.  J.  3.  Hudson 

Materials  Division 

Rensselaer  Polytechnic  Institute 

Troy,  New  York  12181  1 

Dr.  John  T.  Yates 
Department  of  Chemistry 
University  of  Pittsburgh 
Pittsburgh,  Pennsylvania  15260  1 

Dr.  Theodore  E.  Madey 

Surface  Chemistry  Section 

Department  of  Commerce 

National  Bureau  of  Standards 

Washington,  D.C.  20234  1 

Dr.  J.  M.  White 
Department  of  Chemistry 
University  of  Texas 

Austin,  Texas  78712  1 

Dr.  Keith  H.  Johnson 
Department  of  Metallurgy  and  Materials 
Science 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139  1 

Dr.  J.  E.  Demuth 
IBM  Cor port ion 

Thornes  J.  Watson  Research  Center 
?.0.  Box  218 

Yorktown  Heights,  New  York  10598  1 


No. 

Copies 

Dr.  C.  P.  Flynn 
Department  of  Physics 
University  of  Illinois 

Urbana,  Illinois  61801  1 

Dr.  W.  Kohn 
Department  of  Physics 
University  of  California 
(San  Diego) 

LaJolla,  California  92037  1 

Dr.  R.  L.  °erk 
Director,  Center  of 
Materials  Research 
University  of  Maryland 

College  Park,  Maryland  20742  1 

Dr.  W.  T.  Peria 
Electrical  Engineering 
Department 

University  of  Minnesota 

Minneapolis,  Minnesota  55455  1 

Dr.  Chia-wei  Woo 
Department  of  Physics 
Northwestern  University 

Evanston,  Illinois  60201  1 

Dr.  D.  C.  Mattls 
Polytechnic  Institute  of 
New  York 
333  Jay  Street 

Brooklyn,  New  York  11201  1 

Dr.  Robert  M.  Rexter 

Department  of  Chemistry 

University  of  Minnesota 

M  nneapolis,  Minnesota  55455  1 

Dr.  R  P.  Van  Duyne 

Chemistry  Department 

Northwestern  University 

Evanston,  Illinois  60201  1 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  056 


No . 

Copies 


No. 


Copies 


Dr.  S.  Sibener 
Department  of  Chemistry 
James  Franck  Institute 
5640  Ellis  Avenue 

Chicago,  Illinois  60637  1 

Dr.  M.  G.  Lagally 
Department  of  Metallurgical 
and  Mining  Engineering 
University  of  Wisconsin 
Madison,  Wisconsin  53706  1 

Dr.  Robert  Gomer 
Department  of  Chemistry 
James  Franck  Institute 
5640  Ellis  Avenue 

Chicago,  Illinois  60637  1 

Dr.  R.  G.  Wallis 
Department  of  Physics 
University  of  California,  Irvine 
Irvine,  California  92664  1 

Dr.  D.  Ramaker 
Chemistry  Department 
George  Washington  University 
Washington,  D.C.  20052  1 

Or.  P.  Hansma 
Chemistry  Department 
University  of  California, 

Santa  Barbara 

Santa  Barbara,  California  93106  1 

Dr.  P.  Hendra 

Chemistry  Department 

Southampton  University 

England  S09JNH  1 

Professor  P.  Skell 
Chemistry  Department 
Pennsylvania  State  University 
University  Park,  Pennsylvania  16802  1 

Dr.  J.  C.  Hemminger 
Chemistry  Department 
University  of  California,  Irvine 
Irvine,  California  92717  1 


Dr.  Martin  Fleischmann 
Department  of  Chemistry 
Southampton  University 
Southampton  S09  5NH 

Hampshire,  England  1 

Dr.  J.  Osteryoung 
Chemistry  Department 
State  University  of  New 
York  at  Buffalo 

Buffalo,  New  York  14214  1 

Dr.  G.  Rubloff 
I.B.M. 

Thomas  J.  Watson  Research  Center 
P.  0.  Box  218 

Yorktown  Heights,  New  York  10598  1 

Dr.  J.  A.  Gardner 

Department  of  Physics 

Oregon  State  University 

Corvallis,  Oregon  97331  1 

Dr.  G.  D.  Stein 

Mechanical  Engineering  Department 

Northwestern  University 

Evanston,  Illinois  60201  1 

Dr.  K.  G.  Spears 

Chemistry  Department 

Northwestern  Uriversity 

Evanston,  Illinois  60201  1 

Dr.  R.  W.  Plummer 
University  of  Pennsylvania 
Department  of  Physics 
Philadelphia,  Pennsylvania  19104  1 

Dr.  E.  Yeager 

Department  of  Chemistry 

Case  Western  Reserve  University 

Cleveland,  Ohio  41106  2 

Professor  D.  Hercules 

University  of  Pittsburgh 

Chemistry  Department 

Pittsburgh,  Pennsylvania  15260  1 


2 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  056 


No. 


Copies 


Professor  N.  Winograd 

The  Pennsylvania  State  University 

Department  of  Chemistry 

University  Park,  Pennsylvania  16802  1 

Profetior-I^  F.  George 
The  University^&f^Rochester 
Chemistry  Department , 

Rochester,  New  York  14627  1 

Professor  Dudley  R.  Herschbach 

Harvard  College 

Office  for  Research  Contracts 

1350  Massachusetts  Avenue 

Cambridge,  Massachusetts  02138  1 

Professor  Horia  Metiu 
University  of  California, 

Santa  Barbara 
Chemistry  Department 
Santa  Barbara,  California  93106  1 

Professor  A.  Steckl 
Rensselaer  Polytechnic  Institute 
Department  of  Electrical  and 
Systems  Engineering 
Integrated  Circuits  Laboratories 
Troy,  New  York  12181  1 

Professor  R.  D.  Archer 
University  of  Massachusetts 
Chemistry  Department 
tonherst,  Massachusetts  01003  1 

Dr.  A.  C.  Pastor 

Hughes  Research  Laboratories 

3011  Malibu  Canyon  Road 

Malibu,  California  90265  1 


3 


